Constrained optimization of sequentially generated entangled multiqubit states 
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We demonstrate how the matrix-product state formalism provides a flexible structure to solve the 
constrained optimization problem associated with the sequential generation of entangled multiqubit 
states under experimental restrictions. We consider a realistic scenario in which an ancillary system 
with a limited number of levels performs restricted sequential interactions with qubits in a row. 
The proposed method relies on a suitable local optimization procedure, yielding an efficient recipe 
for the realistic and approximate sequential generation of any entangled multiqubit state. We give 
paradigmatic examples that may be of interest for theoretical and experimental developments. 

PACS numbers: 03.67.Bg, 02.70.-c, 42.50.Dv, 71.27.+a 



I. INTRODUCTION 

Entangled multiqubit states are of central importance 
in the fields of quantum computation and quantum com- 
munication [it], and have been the subject of intensive 
theoretical and experimental investigations. As pointed 
out by Schon et al. 0, S| , the classes of all sequentially 
generated multiqubit states, assisted by an itinerant an- 
cilla, are exactly given by the hierarchy of matrix-product 
states (MPSs) y. In this context, the required num- 
ber of ancilla levels is determined by the dimension of 
the MPS canonical representation of the target multi- 
qubit state. Matrix-product states play an important 
role in the context of strongly correlated systems [H| 
and describe the approximate ground states produced 
by density-matrix renormaHzation group (DMRG) qJ0| 
and Wilson's numerical renormalization group |9||. 
Paradigmatic multiqubit states, such as Greenberger- 
Horne-ZeiHnger (GHZ) W [lH, and cluster |l3| 

states, can be described by low-dimensional MPS and 
are considered valuable resources for quantum informa- 
tion and communication tasks. 

The generation of multiqubit entangled states via a 
single global unitary operation acting on initially decou- 
pled qubits is in general a difficult problem. From this 
point of view, several theoretical and experimental ef- 
forts have been oriented towards the sequential genera- 
tion of paradigmatic entangled multipartite states. As 
a matter of fact, a number of sequential and global ap- 
proaches have been implemented in different ph;Y^ical sys - 
tems to produce specifically GHZ Hill, [iB,[ii,[i3|, 
and cluster [3| states. In order to generate sequentially 
any multiqubit state, a wide range of ancilla levels and 
ancilla-qubit operations are necessary In this sense, 
two important theoretical and experimental questions ap- 
pear naturally: will the sequential generation of a desired 
multiqubit state still be feasible under given restricted ex- 



perimental conditions? And if the answer is no, can we 
design an efficient protocol that tells us the best possi- 
ble approximation to the sequential generation of such a 
state? In this paper, we answer satisfactorily both ques- 
tions. We demonstrate how the MPS formaHsm allows 
us to exploit linear algebraic tools to study this relevant 
constrained optimization problem [l9| . 



II. RESTRICTIONS ON THE NUMBER OF 
ANCILLA LEVELS 



It is known that any n-qubit state can be writ- 
ten canonically as an MPS with minimal dimension D{< 
2") It was also shown that such a state can be built 
sequentially with a D-dimensional ancillajif we have ac- 
cess to arbitrary ancilla-qubit unitaries [2|. In the se- 
quential generation of states, an ancillary system A (e.g., 
a Z3-level atom) couples sequentially to an initially decou- 
pled qubit chain jV'/) = IV'i"') ® • ■ ■ ® IV"/^') (e.g., cavity 
photonic qubits that leak out after interacting with an 
atom). Assuming that in the last step the ancilla de- 
couples unitarily from the multiqubit system, we are left 
with the n-qubit state Q 
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.V,'Mipi)\ln,...,ll) , (1) 



an MPS of bond dimension dim(|?/;}) = D, where the 
{D X _D)-dimensional matrix Vj^^ represents the ancilla- 
qubit operation at step k of the sequential generation 
(with isometry condition Yl\k=o V^'^^V'^'' — 1), with \ipi) 
and \'^f) being the initial and the final ancilla states, re- 
spectively. Hence, a relevant experimental question may 
be raised: how well can we represent a given multiqubit 
state I V') if only an ancilla with a smaller number of levels, 
D' < D, is available? More formally: given a state 
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with a canonical MPS representation of bond dimension 
D, what is the optimal MPS of lower bond dimen- 
sion D' < D that minimizes their distance? We want to 
estimate 



min II IV') - 1^) r 

dim{\4>))=D' <D 



(2) 



We propose two techniques to perform the MPS approx- 
imation above, both exploiting a suitably designed local 
optimization of the V matrices in Eq. (U). In the first ap- 
proach, we make use of a corollary of the singular value 
decomposition (SVD) theorem from linear algebra to per- 
form a local optimization procedure which may be called 
"MPS compression," in analogy to the image compression 
tech niq ue already used in computer science and engineer- 
ing [20(1 . Let the SVD of matrix A with rank(A) = r be 
given by ^ = J2i=i o'iUivj. Then, the best possible lower- 
rank approximation to A that minimizes the Frobenius- 
norm distance min^^^j^j^^^)^^,^^ \\A - A\\p is given hy A ^ 

^[^-^ (TiMiuJ [2ll. [2^. This suggests a truncation scheme 
in which one keeps only the r' largest singular values of 
A to form the optimal lower-rank matrix A. We exploit 
now this property, valid for a single matrix, and apply 
the outlined truncation to each matrix Vj^^p k — 1, 
in Eq. ([l]), yielding an MPS of lower bond dimension 
D' = D — {r — r'). This method offers a good solution 
for matrices with well-decaying singular-value spectrum. 

In the second approach a DMRG-inspired varia- 
tional optimization of V matrices [l^l, we seek the best 
possible approximation to \^) in the space of all MPS 
IV") of the form |[T]) (with V ^ V) with bond dimension 
D' < D,hy solving the minimization problem of Eq. ^ 
under the constant-norm condition (V'lV') = 1j which is 
implemented using a Lagrange multiplier A. Varying 
Eq. |[2|) with respect to the matrices defining |V') leads 
to a set of equations, one for each iu, of the form 
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(1 + A) (^|^)-2Re((^|^)) 



= 0, 



(3) 



which determines the optimal V matrices of the desired 
state |-0). These equations can be solved very efficiently 
using a "sweeping procedure" in which one fixes all but 
the fcth ^-matrix and solves the corresponding Eq. ^ for 
the matrix V^^'j . Then one moves on to the neighboring 
site and, in this fashion, sweeps back and forth through 
the chain until the convergence is reached. 

Figure [1] illustrates the two optimization schemes out- 
lined above for two different states, both with D = 16, 
namely (i) the ground state of the XXZ Heisenberg 
Hamiltonian and (ii) a randomly chosen MPS. For (i), 
which has a well-decaying singular-value spectrum, the 
ancilla dimension can be effectively reduced from 16 to 
6. Since variational optimization allows for the feedback 
of information by several sweeps, it generally performs 
better than MPS compression. 
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Figure 1: (Color online) Comparison of the variational op- 
timization approach (solid lines) with the MPS compression 
technique (dotted lines). We consider the ground state of 
the XXZ Heisenberg Hamiltonian (circles) and a randomly 
initialized MPS (triangles), indicating how well these MPS 
with bond dimension D can be approximated with those of 
dimension D' < D. 



III. RESTRICTIONS ON THE SOURCE-QUBIT 
INTERACTIONS 



Every open-boundary MPS of form |(T|) (with V ^ A) 
with arbitrary A matrices, not necessarily isometrics, can 
be cast into a canonical MPS representation with mini- 
mal dimension D [23|. Such states, as mentioned above 
Eq. (dJ, can be generated sequentially 0, such that the 
ancilla decouples unitarily in the last step. We note that 
the sequential generating isometrics can be constructed 
explicitly by successive SVD of the A matrices and ex- 
ploiting the gauge freedom of the matrix-product states 
as outlined in Refs. 0, 01 ■ This is a general recipe for 
the sequential generation of an arbitrary entangled mul- 
tiqubit state if the required ancilla dimension D and 
ancilla-qubit unitaries are available. However, in gen- 
eral, a given physical setup may not have access to some 
of the required local ancilla-qubit unitaries. Given such 
a limitation, we face an interesting constrained optimiza- 
tion problem: which is the sequential protocol by which 
a given multiqubit "target" state can be approximately 
generated with a maximal fidelity? 

To address this problem, let us begin by considering 
the general unrestricted case: the unitary time evolution 
of the joint system ancilla-qubit at step k of the sequen- 
tial generation may be described by a general unitary 



AB _ 

[k] 



w I where 



is a general bipartite Hamiltonian that couples the 
ancilla with the fcth qubit. The latter can be written as 



ttAB 

^[k] 



J2L,,B=ohfl,^<JjA<»'Jj, where hf^^^ are real- 



3AJB=0 JaJB 1-^ 

valued coupling constants and ci, (T2, are the usual 
Pauli a matrices, with ctq = / as the identity matrix. For 
the sake of simplicity, we have considered the case D = 2, 
but similar generators can be found for D > 2. 

Now, suppose that only a restricted set of unitaries 
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Figure 2: The contraction pattern used to calculate the cost 
function in Eq. ^ including the local ancilla operations lA'^ 
and local qubit operations U'^ . The initial states of the qubits 
are denoted by 
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are available. As an illustrative case, let the entan- 
gling Hamiltonian have the restricted form of the XY 
model liil 

^[fcf = hf{a^®cy^+a^®cy^), (4) 

containing a single nonzero contribution = hf'j = 

/122' • Given an arbitrary MPS of the form of Eq. |(T]) (with 
V ^ A) with arbitrary A matrices and the restricted 
Hamiltonian of Eq. |[4]), the aim is to find the optimal re- 
stricted unitary operations Uj^^ = e"*^!*"! that when 
applied sequentially to an arbitrary initial state of the 
joint system |$/) = \(pi) yield a state of the form 

l^) = t>[^f ■••C>[^f?^[ffl*/), (5) 

which is "closest" to the target state of the form |(/?f} (X" 
V'}, where \ipf) is arbitrary. Note that the action of each 
restricted unitary on initial state of qubit, U^^\ip^l^^), 
produces a restricted isometry of the form 

with the definition f^^^ = J^j, K' fijklijjf^) for the re- 
sulting isometry Vik]^ ■ In the ideal case, when the fidelity 
reaches unity, the ancilla can be set to decouple unitar- 
ily in the last step. However, this will not be the case 
in general when the allowed ancilla-qubit unitaries are 
restricted. Thus, the optimization problem reads 

min II I*) - \ipF) (g> \^) f , (7) 

involving a multivariable cost function in \ipF) and 
{h'^i \ . . . , h^i^}, with h^i^ — h^ih, as the variational pa- 
rameters, which can be solved in an iterative procedure. 
We start by picking a particular unitary, say , and 
minimizing the cost function in Eq. ([7]), varying over 

h[ ' , and regarding coupHngs of all the other unitaries 
as fixed. Then we move on to the neighboring unitary 



Figure 3: (Color online) The deviation of the fidelity 1 — J- = 
1 — IK^IV')!! as a function of the number n of qubits for the 
W state with D = 2 when optimizing the couplings ^jyug 
and the local ancilla unitaries Ua, with initial qubit states all 
equal, \tlj\''^) = |0). The inset shows the case where only the 
couplings hjj^j^ are being optimized. 



and optimize its coupling. When all unitaries have been 
optimized locally, we sweep back again and so forth until 
convergence. Each iteration of the local optimization pro- 
cedure requires the calculation of the overlap of the states 
in the cost function of Eq. ([7]) , which can be straightfor- 
wardly calculated in MPS representation as illustrated 
in Fig. [2] (with and set to 1 there) . Varying over 
the vector \ipp) and using the resulting optimal one, the 
cost function simplifies to 2(1 — || (vE'jV') ||), suggesting the 
definition of the fidelity of the procedure as, T = 

For the restricted entangling Hamiltonian of Eq. 
the variational space is so small (only one parameter at 
each step) , that the variational optimization procedure in 
general does not result in much overlap with the target 
state 1-0), as illustrated in the inset of Fig.[3]using the fa- 
miliar \W)n state as target. However, J- can be improved 
by enlarging the variational space. For example, consider 
[/-4i3 jjj Q being replaced with restricted unitaries of 

the form U^^^f , where U^^ = e"*^w*/'^ are arbitrary 
local ancilla unitaries of dimension DxD. This optimiza- 
tion problem can be treated in the same manner as the 
one described in Eq. except that before optimizing 
each C/-^^, we will also vary over the ancilla operation W-^. 
In this way, we are able to produce the |Vl^)„-state with 
almost perfect fidelity (e.g., 1 — w 10~^ for n = 4) as 
illustrated in Fig. [31 In both cases, the smaller the num- 
ber of qubits n, the larger the fidelity, which is a purely 
numerical issue due to the local optimization. Models 
requiring the entangling Hamiltonian of the XX Z form 

h^i^ (cTi ® (Ti -I- (72 (Xi (72) -I- h^2^a3 (g) 0-3, can be simulated in 
a similar manner. 

As a test of the proposed protocols, we appHed our 
variational prescription to the sequential generation of 
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Table I: Comparing the optimal couplings of our simulation hf"^ to those used for experimental realization of W state 
in Ref. [l^l for n — 5. 

Site index (fc) 1 2 3 4 5 



[(/if V^r''*) - 1] X lO'^ 36.50 0.72 8.64 0.62 0.59 



W states in an ion chain. Following closely the recent 
experiment of Ref. , we targeted a W state with the 
entangling Hamiltonian of the form /ii((T+ ® (7+ + ct^ (g) 
a~), with cr+ and a~ being the usual raising and lowering 
Pauli operators, respectively, and the initial state \ipi) ~ 
1)|0) . . . |0) used in experiment. The optimal couplings 
hf"^ of the resulting converged variational MPS \^) (with 
1—T « 10"^ for n = 5) turned out to agree very well with 
the two-qubit rotations h^^^^ used for the experiment of 
Ref. fUl, as illustrated in Table H 

As the main result of this paper, we have found strong 
numerical evidence that an arbitrary MPS with D — 2 
can be generated sequentially if the single-parameter re- 
stricted unitaries If-^'^ in Eq. |[5|) [based on Eq. ^] are 
augmented by arbitrary local unitaries for both ancilla 
and qubit spaces. The combined unitary employed was 
^[fc]^[fc] ^[fef ^[fcf ' where U^^^ = e-'"w^/^ are arbitrary lo- 
cal qubit unitaries (see Fig. [2]). We have considered, for 
this purpose, the generation of 100 randomly chosen MPS 
and have found that 1 — J- remains below 6 x 10^^'^ up to 
n = 5. Note that the combined action of these unitaries 
includes (at most) 11 real independent parameters, which 
in practice can be reduced to ten, since varying a global 
phase has no effect. In contrast, the unrestricted uni- 
taries U-^'^ involve 16 real independent parameters. Thus 
sequential generation of an arbitrary MPS with D = 2, 
can be achieved more economically than previously real- 
ized: a sufficient condition is the availability of the set of 
restricted two-qubit isometries specified above, instead 
of the availability of arbitrary two-qubit unitaries Pi- 

We may then wonder whether some fixed parameter- 
free two-qubit isometries can act as universal set for gen- 
eration of arbitrary entangled states. The problem we 
propose, which is the natural one in the sequential gen- 
eration of multi-qubit states, is the following: give a 
minimal set S of two-qubit unitaries such that one can 
generate an arbitrary isometry with a single unitary of 
the set S, together with arbitrary one-qubit unitaries. 
Note that, we already showed numerically that S can 
be given by the single-parameter interactions of the XY 
type, whereas we now wonder whether this can be re- 
alized by a minimal set of fixed canonical gates. Note 
that since the paradigm is completely different (a single 
use of the entangling unitary and isometries instead of 
unitaries), the results concerning universal sets of gates 
for quantum computing do not play a role for our proto- 
col. We have found numerically, for example, that some 
parameter-free fixed two-qubit gates [such as controlled 



NOT (CNOT)] plus three local unitaries are not isomet- 
rically universal, as they are not capable of generating 
an arbitrary state with J- — 1. The search for such two- 
qubit gates, if any, remains open. 

Recently, a lot of effort has been devoted to find min- 
imal sets of one-qubit and two-qubit gates, and the min- 
imal number of applications, to generate arbitrary two- 
qubit unitaries [2Q[. The existence of these universal sets 
is of central relevance in quantum computing. The above 
results suggest consideration of a class of problems in- 
volving a different paradigm: which are the universal sets 
of one-qubit and two-qubit gates that can generate arbi- 
trary two-qubit isometries? What is the minimal num- 
ber of applications and how does this compare to the 
quantum computing case? For the case of two-qubit uni- 
taries, a universal gate set (in the usual quantum com- 
puting sense) is clearly sufficient, but not necessary. This 
results, for example, from counting the number of inde- 
pendent parameters for an arbitrary two-qubit unitary, 
clearly larger than in the case of an arbitrary two-qubit 
isometry. The aim will be then to find the exact de- 
composition of an arbitrary isometry into a minimal ap- 
plications of unitaries as computational primitives. The 
general solution associated with this paradigm remains 
open. 

Finally, we also want to point out that our scheme by 
construction can be clearly viewed also within the general 
framework of optimal control theory (27l . [28| . 



IV. CONCLUSIONS 

In conclusion, we have developed protocols for an effi- 
cient sequential generation of entangled multiqubit states 
under realistic experimental constraints. We stress that 
the proposed optimization methods are of wide applica- 
bility and will be of importance for any sequential physi- 
cal setup. In particular, we can mention photonic qubits, 
atoms, ions, superconducting qubits, or quantum dots. 
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